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Abstract 

We used microsatellite markers to assess the population genetic structure of the 
scribbled rabbitfish Siganus spinus in the western Pacific. This species is a cultur- 
ally important food fish in the Mariana Archipelago and subject to high fishing 
pressure. Our primary hypothesis was to test whether the individuals resident in 
the southern Mariana Island chain were genetically distinct and hence should be 
managed as discrete stocks. In addition to spatial sampling of adults, newly-settled 
individuals were sampled on Guam over four recruitment events to assess the tem- 
poral stability of the observed spatial patterns, and evidence of self- recruitment. 
We found significant genetic structure in S. spinus across the western Pacific, with 
Bayesian analyses revealing three genetically distinct clusters: the southern Mariana 
Islands, east Micronesia, and the west Pacific; with the southern Mariana Islands 
being more strongly differentiated from the rest of the region. Analyses of temporal 
samples from Guam indicated the southern Mariana cluster was stable over time, 
with no genetic differentiation between adults versus recruits, or between samples 
collected across four separate recruitment events spanning 1 1 months. Subsequent 
assignment tests indicated seven recruits had self- recruited from within the South- 
ern Mariana Islands population. Our results confirm the relative isolation of the 
southern Mariana Islands population and highlight how local processes can act 
to isolate populations that, by virtue of their broad- scale distribution, have been 
subject to traditionally high gene flows. Our results add to a growing consensus 
that self-recruitment is a highly significant influence on the population dynamics 
of tropical reef fish. 



Introduction 

For most marine organisms, a pelagic larval stage provides 
the primary mechanism for dispersal among spatially frag- 
mented habitat patches (Kritzer and Sale 2004; Co wen and 
Sponaugle 2009). The degree to which larvae disperse and 
populations are connected has a profound influence on the 
population dynamics of a species (Hixon et al. 2002; Gaines 
et al. 2007), with the stability and resilience of populations 
dependent upon a constant supply of larvae, either locally 
or externally sourced (Warner and Cowen 2002). Despite the 
obvious need, quantification of connectivity has proven dif- 
ficult. This is principally due to the small size of larvae, their 



patchy distribution, and high rates of larval mortality (Leis 
1991), all of which has severely hindered attempts at in situ 
studies. 

Because the pelagic larval duration of many fish larvae 
can be more than one month (Victor and Wellington 2000), 
coral reef fish populations were assumed to be demograph- 
ically open with a large dispersive potential (Roughgarden 
et al. 1985; Caley et al. 1996). In addition, studies on 
settlement- stage larval fish behavior have shown that many 
species are strong swimmers (Stobutzki and Bellwood 1997; 
Fisher 2005) and can respond to a variety of potential settle- 
ment cues, both olfactory (Atema et al. 2002; Lecchini et al. 
2005) and auditory (Leis et al. 2002; Simpson et al. 2005). 
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The combination of these traits was thought to favor maxi- 
mization of an individual's dispersal distance. However, over 
the past 10 years studies have shown that these same traits 
can equally be used to actively reduce dispersal and entrain 
larvae close to natal reefs (Gerlach et al. 2007). 

Traditionally, the most widely used approach for estimat- 
ing connectivity has been population genetics (see reviews by 
Hellberg 2007; Jones et al. 2009). The analysis of genetic vari- 
ation among spatially isolated populations allows for an in- 
direct assessment of connectivity (Neigel 1997). However, in 
many marine organisms gene flow is high over evolutionary 
timescales, with only a few successful migrants per generation 
needed to produce genetic homogeneity among populations 
(Slatkin 1993). This results in connectivity estimates that re- 
flect historical processes over multiple generations (Hellberg 
2007). Additionally, population-level estimates of connectiv- 
ity require several simplif)^ing assumptions, which may be 
hard to satisfy in natural biological systems (Whitlock and 
McCauley 1999). Such limitations have recently been over- 
come with the development of new, more powerful anal- 
yses using highly variable micro satellite markers that have 
shifted the focus of analysis from populations to individuals 
(Ballouxand Lugon-Moulin 2002; Pearse and Crandall 2004). 
In addition to using more variable markers, advances in statis- 
tical methodologies using maximum likelihood and Bayesian 
inference enables sorting of genetically similar individuals 
into discrete populations, and also enables individuals to be 
assigned to their population of origin (Man el et al. 2005), 
thus enumerating rates of contemporary connectivity (e.g.. 
Underwood et al. 2007; Saenz-Agudelo et al. 2009). Direct es- 
timates can also be obtained using parentage analysis (Planes 
et al. 2009) , which has been shown to provide estimates of self- 
recruitment identical to tagging studies (Jones et al. 2005). 
While parentage analysis was once restricted to fish species 
in which all prospective parents could be sampled (see e.g., 
Jones et al. 2005; Planes et al. 2009), novel Bayesian parentage 
methods now make it possible to assess parentage when only 
a small proportion of potential parents are sampled (Christie 
2010). 

Utilizing micro satellite markers and advances in statistical 
capabilities, research efforts in recent years have seen our 
understanding of population connectivity in reef fish increase 
dramatically. Studies have documented high levels of self- 
recruitment in species with large potential dispersive capacity 
(Jones et al. 1999, 2005; Almany et al. 2007). No populations, 
however, have been found to be completely closed (Jones 
et al. 2009). Larvae have been shown to disperse up to 35 km 
away from natal sites in populations with high local retention 
rates (Planes et al. 2009), and extensive gene flow between 
populations is commonly observed (van Herwerden et al. 
2003; Purcell et al. 2006; Home et al. 2008; Christie et al. 
2010). Consequently, coral reef fish populations cannot be 
categorized as merely open or closed (Co wen et al. 2000; 



Mora and Sale 2002) with many populations appearing to 
exhibit both self-recruitment and long-range dispersal, the 
ratio of which may vary dynamically with location and time 
(Cowen and Sponaugle 2009). 

To date, the focus of most connectivity studies has been 
quantification of spatial patterns with little knowledge of the 
temporal stability of demo graphically relevant connectivity 
(but see Selkoe et al. 2006; Jones et al. 2010). However, dis- 
crete cohorts can possess unique genetic signatures (Planes 
and Lenfant 2002) and examination of the variability of these 
genetic signatures over time can lead to insights into the na- 
ture of larval dispersal. One process that may dramatically 
affect temporal gene flow is sweepstakes reproductive suc- 
cess, where only a small proportion of reproductive adults 
are responsible for the majority of the recruitment within a 
population (Li and Hedgecock 1998; Hedgecock et al. 2007). 
This reduced genetic contribution to subsequent recruitment 
leads to testable hypotheses regarding genetic diversity and 
relatedness within and among recruits and adult populations 
(Hedgecock 2010). Analysis of multiple dispersal events are 
needed to resolve and understand the complex larval disper- 
sal patterns of coral reef fish, which can ultimately be used to 
increase our understanding of population dynamics and lead 
to better resource management (Selkoe et al. 2006; Planes 
et al. 2009; Christie et al. 2010). 

In this study, we investigate the population structure of 
the scribbled rabbitfish {Siganus spinus) in the western Pa- 
cific, with a particular focus on the islands of the southern 
Mariana archipelago (Fig. 1). This herbivorous species is a 
highly fecund, fast- growing, common inhabitant of shallow 
coral reefs throughout the Indo-Pacific region (Woodland 
1990). On Guam, the largest and most southerly island in the 
Mariana chain, settlement-stage S. spinus recruit to shallow 
reef flats in large numbers, usually twice a year (Tsuda and 
Bryan 1973), at the end of a 32-day pelagic phase (Chirichetti 
1996). These recruits settle at a relatively large size, 




Figure 1. An adult scribbled rabbitfish (Siganus spinus) pictured at night 
in Piti Marine Preserve, Guam. Photo by Mark Priest. 
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approximately 43-mm fork length (FL) (Kami and Ikehara 
1976), and have strong swimming capabilities (Fisher 2005). 
Newly settled recruits are considered a cultural delicacy and 
are harvested en masse by local residents. Adult S. spinus 
are also a favored food fish heavily targeted by spear and 
net fishermen. Recent declines in Guam's reef fish stocks 
(Newton et al. 2007; Zeller et al. 2007) and uncertainties 
over population size and stock structure have seen Guam's 
rabbitfish populations recognized as a "species group of 
greatest conservation need" by the local fisheries agency 
(Bassler and Aguon 2006). Should self-recruitment be a 
significant pathway for population replenishment, any de- 
crease in spawner biomass as a result of overharvesting could 
have considerable effects on subsequent recruitment events 
and the long-term sustainability of the fishery (Man et al. 
1995). 

Here, we use six microsatellite markers to investigate (1) 
the spatial genetic population structure of S. spinus at regional 
(western Pacific) to local (southern Mariana islands) scales. 
Prevailing large-scale current patterns indicate an east-west 
flow between islands in this region and analysis of gene flow 
pathways will identify any genetic discontinuities and pos- 
sible barriers to dispersal. The lack of any obvious physical 
barriers would implicate oceanographic conditions in such a 
case. Sampling from multiple locations around Guam pro- 
vides insights into whether significant structure is evident at 
the island scale. Such information is important in the con- 
text of formulating effective conservation management plans 



for this species; (2) in conjunction with the local-scale spa- 
tial sampling we also sampled genetic variability within and 
between multiple cohorts at several sites on Guam, to assess 
the temporal stability of gene flow, degree of self-recruitment, 
and evidence of sweepstakes reproductive success. Many con- 
nectivity studies on coral reef fish that have included temporal 
sampling have focused on small, long-lived site-attached fish 
(see e.g., Hepburn et al. 2009; Planes et al. 2009; Christie et al. 
2010; Jones et al. 2010), whereas this study provides one of 
the first assessments of a more mobile reef fish species. 

Methods 

Sampling regime 

To assess spatial genetic structure at regional scales (200-5000 
km), tissue samples from 971 S. spinus individuals were col- 
lected from nine islands across the western Pacific (Fig. 2; 
Table 1) between August 2007 and November 2009. Addi- 
tionally, where possible, individuals were collected from at 
least two discrete sites at each island to test for population 
structure at local (within island, 1-60 km) spatial scales. In- 
dividuals were sampled across the whole size range of the 
species (>60 mm FL) to capture as much genetic diversity 
of the population as possible, and to minimize the influ- 
ence of any single, potentially genetically atypical, large co- 
hort. To assess temporal genetic structure on Guam, 331 
newly settled S. spinus recruits (<50 mm FL) were collected 
from four discrete recruitment events: November 2007, July 
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Figure 2. Sampling locations for Siganus spinus across the western Pacific. For detailed locations of sites and subsites, see Table 1 . Prevailing regional 
surface current direction is denoted by arrows with speed represented in color. Current data are long-term mean (1993-2009) obtained from the 
NOAA ocean surface current analyses — real time website (http://www.oscar.noaa.gov). 
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Table 1. Summary of Siganus spinas sampling locations, site codes, sample sizes (A/), and genetic diversity measures over all loci for the large-scale 
spatial analyses. Observations include mean number of alleles (A/a), number of private alleles (A/ra), allelic richness {AR), observed (Ho) and expected 

(He) heterozygosity, and mean inbreeding coefficient (Fis). 
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2008, September 2008, and October 2008 from nine locations 
(Fig. 3; Table 2). The predictable timing of recruitment 
(±3 days over the last quarter moon; Kami and Ikehara 1976) 
allowed separate cohorts of recruits to be clearly identified as 
belonging to discrete recruitment events. However, the vari- 
able nature of recruitment meant individuals could not be 
sampled from every site for each recruitment event. 

All fish were collected from shallow reef- flats using hand 
spears, with the exception of the Philippine samples, which 
were obtained from the Dumaguete fish market. Individuals 
were kept in ice slurry before being processed. Tissue samples, 
taken as fin clips, were preserved in 95% ethanol (EtOH) and 
stored at room temperature while awaiting further analysis. 

DNA extraction and laboratory analysis 

Genomic DNA was extracted using DNeasy Blood and Tissue 
kits (QIAGEN, Valencia, CA, USA). Samples were genotyped 
at five micro satellite loci originally designed for the S. spinus 
congener S.fuscescens (Ravago-Gotanco et al. 2010). The loci 
used in this study were Sfus-8, Sfus-9, Sfus-56, Sfus-98, and 
Sfus-113. Samples from Guam were also genotyped for the 
additional marker Sfus-21. Polymerase chain reactions (PGR) 
were carried out individually in a 10 /xl reaction volume 
containing: 1 /xl lOx Taq polymerase buffer, 0.2 /xl Taq DNA 
polymerase, 0.2 /xl lOmM dNTPs, 0.5 /xl 25mM MgGl2, 5.1 
/xl ddHaO, 1 /xl of genomic DNA (~ 10-40 ng /xP^), and 
1 /xl each of forward and reverse primers. For loci Sfus-8, 



Sfus-21, Sfus-56, and Sfus-113, primer concentrations were 
5 /xM; for locus Sfus-98, primer concentration was 2.5 /xM. 
Thermocycling conditions for the PGRs consisted of an initial 
step of 80 G for 1 min and a denaturing step of 94 G for 
2 min. Ten cycles of 30 sec at 94 G, 30 sec at 56 G, and 30 sec 
at 72 G were followed by 20 cycles of 30 sec at 94 G, 30 
sec at 54- G, and 30 sec at 72 - G, with a final extension step 
of 4 min at 72 G. Gapillary electrophoresis was performed 
on pooled PGR products using an AB 3730 DNA Analyzer 
(Applied Biosystems, Garlsbad, GA, USA), and scored with 
GeneMapperv. 3.7 (Applied Biosystems, Garlsbad, GA,USA). 
All PGR reactions and genotyping were performed by the 
Australian Genome Research Facility, Brisbane, Australia. 

Genetic analyses 

The mean number of alleles per locus, allelic richness, and 
number of observed alleles were calculated using fstat v. 
2.9.3 (Goudet 200 1 ) . Expected and observed heterozygosities, 
and the number of private alleles per population, were eval- 
uated using GenAlEx v. 6.1 (Peakall and Smouse 2006). De- 
partures from Hardy- Weinberg equilibrium (HWE) within 
all populations by locus, and across all loci were examined 
using GENPOP V. 3.4 (Raymond and Rousset 1995); signif- 
icance levels were adjusted with sequential Bonferroni cor- 
rection for multiple tests for P < 0.05. To assess the cause 
of any deviations from HWE, micro-checker (van Oost- 
erhout et al. 2004) was used to assess the data for null 
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Figure 3. Sampling sites for adult arid recruit collectioris of Siganus 
spinas from Guam. See Table 2 for sample size and time of collections. 
Sampling was not conducted on the northeast of Guam as there is no 
suitable reef flat habitat. 



alleles, stuttering, and large-allele dropout. Loci were tested 
for linkage disequilibrium using genpop (1000 batches, 1000 
iterations) for each locus pair across all populations. One 
hundred and twenty-eight samples were regenotyped to de- 
termine the study-specific error rate and resolve suspected 
null homozygotes. 



Large-scale population structure and 
assignment tests 

Population structure was assessed using standard Fsr genetic 
differentiation measures and Bayesian assignment methods. 
Global and population pairwise Fsr values between all sam- 
pled sites were calculated using FreeNA (Chapuis and Estoup 
2007), with and without correction for null alleles, and as- 
sessed using Fisher's exact tests of significance including se- 
quential Bonferroni correction for multiple tests for P < 0.05 
with FSTAT. Regional genetic structure was assessed at the 
island level using analysis of molecular variance (AMOVA) 
(Excoffier et al. 1992) implemented in arlequin v. 3.1 1 (Ex- 
coffier et al. 2005). Significance of Fsr values was calculated 
using 10,000 Markov Chain Monte Carlo (MCMC) permu- 
tations of alleles across clusters, also using arlequin. Results 
were visualized by performing a principle coordinate anal- 
ysis (PCoA) on a genetic distance matrix constructed using 
pairwise Fsr values. Patterns of isolation-by-distance were 
investigated by plotting genetic distance against geographic 
distance and tested with a mantel test using 999 permuta- 
tions implemented in GenAlEx. Cavalli-Sforza and Edwards' 
(1967) chord distance (Dq) was used as it is insensitive to 
departures from HWE. Geographic distances were calculated 
as the shortest distance by sea between the midpoints of sam- 
pling locations. 

To investigate the presence of subtle genetic structure and 
spatial genetic discontinuities, we used Bayesian assignment 
methods to determine the most likely number of homoge- 
nous genetic clusters of all sampled individuals. First, struc- 
ture V. 2.2.3 (Pritchard et al. 2000; Falush et al. 2007) was 
used to assign individuals to a predetermined number of sub- 
groups {K) based on the likelihood of the individual belong- 
ing to each subgroup, as determined by the individual's allele 
frequencies at each locus. K was set from 1 to 18 and each 
model run was repeated independently 10 times using a run 
length of 100,000 MCMC repetitions and a burn-in period 



Table 2. Summary of Siganus spinas sampling locations, times, and sample sizes for Guam. 
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includes individuals >100 mm FL used in large-scale spatial analyses. 
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of 100,000. The model was run with admixture, correlated 
allele frequencies, and with prior location options (Hubisz 
et al. 2009). The most appropriate K was determined from 
plots of ad hoc posterior probability models of AiC, as ex- 
tracted by STRUCTURE HARVESTER V. 0.56.4 (Evauno et al. 
2005; Earl and von Holdt 2011). Second, the geneland v. 
3.2.1 (Guillot et al. 2005, 2008) package for the R statisti- 
cal software environment (R Development Core Team 2009) 
was employed to assess spatial genetic discontinuities be- 
tween sampling locations, geneland uses a Bayesian cluster- 
ing model similar to structure, but also includes geographic 
sampling location data as a weak prior to produce spatially 
explicit maps of genetically homogenous clusters. Again, 10 
independent runs were used for each Ky with K ranging from 
1 to 18, using 1,000,000 MCMC repetitions with a burn-in 
period of 200,000 and a thinning value of 1000. 

In order to assess the origin of Guam's recruits, we used 
the Bayesian assignment method of Rannala and Mountain 
(1997) implemented in geneclass2 (Piry et al. 2004) using 
the clusters detected by structure as source populations. 
Recruits were assigned to a population if the probability of 
assignment was greater than 0.05 to only one source pop- 
ulation, as determined using the simulation algorithm of 
Paetkau et al. (2004) (10,000 simulations). 

Fine-scale spatial structure and temporal 
comparisons 

To assess fine-scale spatial and temporal patterns of genetic 
differentiation on Guam, AMOVA was used in arlequin 
to partition the genetic variance between location and life- 
history stage, that is, adults and recruits. This evaluates ge- 
netic differences based on allele frequency variation between 
all samples. We also constructed a pairwise matrix of related- 
ness values between all adult and recruit samples from Guam 
using Queller and Goodnight s relatedness metric (Queller 
and Goodnight 1989), Rxy^ implemented in GenAlEx. This 
metric describes the relationship of pairs of individuals based 
on the number of shared alleles standardized by allelic state. 
Adults were classified as individuals >100 mm FL, as this is 
the size at first maturity (unpubl. data). PCoA was used to vi- 
sualize the results and nonparametric multivariate analysis of 
variance (MANOVA) (Anderson 2001; McArdle and Ander- 
son 2001) was used to test whether individuals within sample 
groups shared more alleles with each other than individuals 
of other groups. Significance was assessed using 999 permuta- 
tions, implemented in primer v. 6 (Clarke and Gorley 2006). 
Tests were performed on all adult and recruit samples for dif- 
ferences between sample site and life history, and a subset of 
data, containing only the four sites at which both adult and 
recruit samples were obtained, was reanalyzed. All recruit 
samples were also tested for differences between recruitment 
events using nonparametric MANOVA. Differences in het- 



erozygosity between adult and recruit samples within all sites 
and recruitment events were assessed in est at using 10,000 
permutations. One would expect increased relatedness (i.e., 
greater shared alleles) within cohorts and decreased related- 
ness between cohorts if sweepstakes reproduction is evident, 
as a result of the high variance of an individuals' chance of 
successfully reproducing. Similarly, one would also expect to 
see decreased heterozygosity in recruit populations compared 
to adult populations (Hedgecock 2010). 

Finally, all recruit samples were added to the region-wide 
dataset and PCoA was performed on pairwise Fst values. 
Samples were grouped into the three clusters determined by 
the region-wide analyses, including Guam recruits as a fourth 
cluster, and pairwise comparisons were made between all 
clusters using nonparametric MANOVA in primer. 

Results 

Genetic diversity estimates were generally low for all sam- 
ple groups. Mean number of alleles per locus within sample 
groups varied from 4.6 to 7.8. The total number of alleles 
at each locus was highly variable, ranging from 23 {Sfus- 
8) to 2 {Sfus-56 and Sfus-9). Both Sfus-56 and Sfus-9 were 
monomorphic within some sample groups. Allelic richness 
ranged from 2.61 to 4.56 and was generally lower in the Guam 
recruit samples than the region-wide samples. A total of three 
private alleles were present in the dataset (Table 1). Exact tests 
for linkage disequilibrium across all sample groups found no 
loci pair was in linkage disequilibrium and that loci could be 
treated as independent. Replicate genotyping of 128 sampled 
showed high levels of repeatability (>96% agreement). 

Significant departures from HWE were detected at locus 
Sfus-8 in 36 of the 40 sample groups (Tables Al and A2). Of 
the other five loci, none exhibited more than three depar- 
tures from HWE across all 40 sample groups. Results from 
the program micro-checker suggest that deviations from 
HWE in Sfus-S are likely due to the presence of null alle- 
les, which resulted in heterozygote deficiencies. The effect 
of null alleles on population differentiation estimates was 
assessed by analyzing the data with FreeNA (Chapuis and 
Estoup 2007), with and without correction for null alleles. 
Global Fsr values showed minimal differences between these 
two datasets (noncorrected Fsr = 0.014; null allele corrected 
Fsr = 0.015) and analysis of pairwise Fsr estimates showed 
no differences between either dataset. Bayesian analyses us- 
ing datasets including and omitting Sfus-8 also resulted in 
congruent outcomes, indicating the presence of null alleles 
had no tangible effect on the results. 

Large-scale spatial population structure 

There was considerable genetic differentiation across the 
study region {Fsr = 0.023, P < 0.001) with 41% of pair- 
wise Fsr comparisons significant (Table 3). A significant 
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Figure 4. Relationship between geographic (Km) and genetic distance 
(Dc) of Sigan us spinas among 1 8 sampled sitesacross the western Pacific. 
Points represent comparisons between pairs of sample groups. 
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Figure 5. Principal coordinate analysis of all pairwise Fsj comparisons 

between the 1 8 sites sampled in the western Pacific. Both axes combined 
explain 77.6% of the total variation. See Table 1 for site codes. Sam- 
ples are color coded by island. Dashed rings denote population clusters 
detected by Bayesian analyses. 



isolation-by-distance relationship was detected throughout 
the study region (Fig. 4) (linear regression: = 0.068; 
P = 0.013, Mantel test), with genetic distance increasing with 
greater geographical separation. Moreover, at the within- 
island scale none of the 16 pairwise Fsr estimates was sig- 
nificantly different from zero, indicating little genetic differ- 
entiation within islands. When samples were grouped at the 
island level, AM OVA showed a significant proportion of ge- 
netic variance (1.7%, P < 0.001) was partitioned among 
islands (Table 4A). This was also supported by a PCoA 
of pairwise Fsr values, where most within-island samples 
clustered closer to each other than between-island samples 
(Fig. 5). The PCoA analysis also revealed a significant dif- 



ference between samples from the southern Mariana Islands 
(Guam and Saipan) and the rest of the study area (nonpara- 
metric MANOVA: F = 57.71, P = 0.001), with the southern 
Mariana samples grouping out along PCoA Axis 1. 

Bayesian analyses used to detect genetic clustering re- 
turned results congruent with PCoA. structure revealed 
three genetically distinct clusters in the data (Fig. 6). Fur- 
ther Bayesian -based analyses with geneland produced sim- 
ilar results to structure. Plots of the posterior distribu- 
tion of the estimated number of populations indicated three 
populations were present in the data, with geneland return- 
ing the same clustering pattern as structure (Fig. 7A-C): 



Table 4. Summary of the analysis of molecular variance used to examine spatial and temporal genetic structure of Siganus spinas. 





Source of variation 


S.S. 


df 


% 


P 


Large-scale spatial comparison 












(A) Grouped by island 


Among Islands 


41.91 


8 


1.74 


<0001 




Among sites within islands 


12.22 


9 


0.02 


0.425 




Among individuals within sites 


1272.48 


953 


22.20 


<0001 




Within individuals 


818.50 


971 


76.04 


0001 


Fine-scale spatial and temporal comparisons 










(B) Guam samples grouped by 


Among adults and recruits 


1.35 


1 


0.00 


0.554 


life-history stage 


Among sites within adults and recruits 


30.31 


20 


0.00 


0.498 




Among individuals within sites 


674.50 


445 


17.11 


<0001 




Within individuals 


501.00 


467 


82.89 


<0001 


(C) Recruits grouped recruitment 


Among recruitment events 


4.14 


3 


0.00 


0.784 


event 


Among sites within recruitment events 


20.76 


13 


0.18 


0.368 




Among individuals within sites 


473.89 


314 


18.06 


<0001 




Within individuals 


346.50 


331 


81.76 


<0001 



S.S., sum of squares; df, degrees of freedom; %, amount of variation explained by specific component; P indicates significance (values <0.05 are 
shown in bold). 
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Figure 6. Large-scale spatial population structure detected \n Siganus 
spinas using the structure program. The peak in AK (open circles) cor- 
responds to the number of homogenous populations {K) detected from 
10 independent runs of each model. AK is an ad hoc statistic based 
on the rate of change of log-likelihood as K is increased. Mean LnP(D) 
(closed circles) is the mean of 10 independent runs. 

Guam and Saipan (southern Mariana Islands); Chuuk, Pohn- 
pei, and Majuro (East Micronesia); and the PhiUppines, 
Palau, Yap, and PNG (West Pacific). This separation into 
three groups is also evident in the pairwise Fsr ordination 
(Fig. 5). The southern Mariana Island samples separate from 
all other samples along PCoA Axis 1, which explains 62.2% 
of the variation in the data, whereas the East Micronesia and 
West Pacific groupings separate out along PCoA axis 2, which 
explains 15.4% of the variation. This pattern is consistent with 
all analyses, suggesting the southern Mariana Islands is the 
most genetically differentiated cluster of three populations 
(southern Mariana Islands, East Micronesia, and the West 
Pacific) present in the study region. 

Assignment tests showed none of 331 individual recruits 
sampled from Guam could be excluded from originating 
from the southern Mariana Islands population, however, 
only seven individuals (2%) were also excluded from both 
the East Micronesia and West Pacific populations, and thus, 
could be assigned as originating from the southern Mari- 
ana Islands population (Table 5). The East Micronesia pop- 
ulation was excluded as a population of origin for recruits 
more often than the West Pacific (9% vs. 2%), potentially 
indicating a greater barrier to dispersal between the South- 
ern Mariana Islands and East Micronesia than the West 
Pacific. 

Fine-scale spatiotemporal population 
structure 

Results indicate genetic homogeny across Guam and Saipan 
both spatially and temporally. Comparisons using AMOVA 
did not reveal any genetic differentiation between adults and 
recruits (Table 4B). PCoA and nonparametric MANOVA on 

(O 2012 The Authors. Published by Blackwell Publishing Ltd. 
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Figure 7. Population structure of Siganus spinas in the western Pacific 
detected using the spatially explicit clustering model in geneland. Each 
panel (A-C) represents the posterior probabilities of membership of one 
of the four populations detected in the analysis. White to red indicates 
high to low probability of membership of the specific cluster. 

pairwise Rxy values for adults and recruits also showed no 
significant differences between life-history stages {F = 0.11, 
P = 0.71), or sample sites (F = 0.62, P = 0.74). There was also 
no evidence of sweepstakes reproduction. Permutation tests 
in FSTAT comparing mean heterozygosity and Fsr between 
adult and recruit samples revealed no differences between 
adult and recruit groups {Hq: P = 0.08; Fsr'- P = 0.51). When 
considering only the recruit samples, there was no difference 
between the eight sampling sites {F = 1.31, P = 0.35). There 
was also no evidence of differentiation between cohorts that 
settled at different recruitment events using either AMOVA 
(Table 4C) or pairwise Rxy analyses {F = 2.23, P = 0.14), 
indicating temporal stability of gene flow over the 1 1 -month 
sampling period. 
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Table 5. Probabilities of exclusion for all recruits excluded from one or more of the three populations detected by structure and geneland. Probabilities 
< 0.05 are shown in bold typeface. 

Probability of population membership 



Recruit ID 


Site 


Rec. event 


E. Micronesia 


S. Mariana Is. 


W. Pacific 




A 

AchanQ 


iMoveiiiijei zuu/ 




0 9RA 
U.Z 


0 1 

U. 1 JO 


SS0243 


icii ly 


Mnwpmhpr 7007 

INHJvCIIIUCI Z.UU/ 


0 020 


0 093 


0.1 22 


550224 


Agat 


Mnwpmhpr 7007 


0 006 


0 342 


0.1 03 


550231 


Agat 


Mnvpmhpr 7007 


0 004 


0 300 


0.027 


550905 


Agat 


Qpntpmhpr 700R 


0 007 


0 355 


0.077 




East Agana 


iMoveiiiijei 1 


U.U 1 13 


0 1 
U.^j 1 


0 1 70 
U. 1 / u 


550947 


CdbL MLJal Id 


JcpLclllUci ZUUo 




0 517 


0 085 


550949 


Fact Anana 
CdjL /-\ydl Id 


Qpntpmhpr 700R 


0 036 


0 427 


0.146 


J jU / / d 


r\\/c^ more 
OUVtrl 1 lUl b 


Qon+prnhpr 9008 
JcpLclllUci ZUUo 


\J.\JJl3 


0 1 37 


0 267 


550786 


OUVtrl 1 lUI J 


Qpntpmhpr 700R 

Jfc:|J Ltrl 1 1 Ucl ZUUO 


0 001 


0 084 


0.002 


550802 


OUVfc;i 1 lUl b 


Qpntpmhpr 900P 
JcpLclllUci ZUUO 




0 327 


0 203 


550286 




Mn\/pmhpr 9007 

IMUvCIIIUCI Z.UU/ 


0 042 


0 1 98 


0.292 


550729 


Pago 


Qpn+pmhpr 900P 
jcpLclllUci ZUUO 


\J.\J \ £. 


0 1 38 


U.U 1 o 


SS0734 


Pago 


September 2008 


0011 


0.604 


0.090 


SS0736 


Pago 


September 2008 


0008 


0.118 


0.026 


SS0755 


Pago 


September 2008 


0042 


0.472 


0.239 


SS0528 


Rios 


July 2008 


0046 


0.316 


0.113 


SS0826 


Rios 


September 2008 


0008 


0.385 


0.038 


SS0829 


Rios 


September 2008 


0015 


0.596 


0.081 


SS0840 


Rios 


September 2008 


0002 


0.192 


0.010 


SS0842 


Rios 


September 2008 


0012 


0.154 


0.073 


SS0848 


Rios 


September 2008 


0028 


0.456 


0.138 


SS0865 


Rios 


September 2008 


0008 


0.057 


0.038 


SS0542 


Tanguisson 


July 2008 


0013 


0.477 


0.108 


SS1028 


Tanguisson 


October 2008 


0049 


0.293 


0.248 


E. Micronesia, 


East Micronesia; S. Mariana Is., 


, southern Mariana Islands; W. Pacific, West Pacific. 







Consistent with no evidence of spatial or temporal ge- 
netic structure within Guam, when the recruit samples from 
Guam were added to the region-wide dataset, the PCoA ordi- 
nation indicated clearly that the Guam recruit samples were 
most closely related to the Southern Marianas Island cluster 
(Fig. 8). Moreover, pairwise nonparametric MANOVA com- 
parisons of recruit samples against the three clusters detected 
by Bayesian analyses showed that the Guam recruit samples 
were not significantly different from the southern Mariana 
Island cluster (t = 1.84, P = 0.133), but were significantly 
differentiated from the two other region-wide clusters (West 
Pacific: t = 7.66, P = 0.001 and East Micronesia: t = 5.53, 
P = 0.001). The similarity between Guam recruit samples 
and the Southern Marianas Island cluster in the region-wide 
analysis is suggestive of limited gene flow between the South- 
ern Marianas Islands and the rest of the region. 

Discussion 

We found significant large-scale spatial genetic structure in S. 
spinus populations across the West Pacific, with the Southern 
Mariana Islands strongly separated from the rest of the region. 
The stability of the Southern Mariana Islands population was 



confirmed by temporal sampling of newly settled fish from 
four separate recruitment events over a period of 1 1 months. 
Recruitment cohorts were genetically indistinguishable from 
each other and from the resident adult population on Guam 
and Saipan. These results indicate limited connectivity with 
the rest of the West Pacific and highlight the likely signifi- 
cance of self- recruitment in maintaining the unique genetic 
signature of reef fish from the Southern Mariana Islands. 

Three genetically distinct populations were found across 
the region: the southern Mariana Islands (Guam and Saipan), 
East Micronesia (Chuuk, Pohnpei, and Majuro), and the West 
Pacific (Philippines, Palau, Yap, and PNG). The Southern 
Mariana Islands were strongly separated from the rest of 
the region in all PsT and Bayesian analyses, indicating sig- 
nificant isolation from the rest of the region. The East Mi- 
cronesian and West Pacific populations exhibited less, but 
nevertheless still significant, differentiation from each other. 
Genetic structure between populations was detected when 
using Bayesian analyses with either prior sampling group 
information (structure) or spatially explicit geographical 
coordinates (geneland). These models have been criticized 
for coaxing data into forming cluster groups by weighting 
results according to sample location. However, in the case of 
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Figure 8. Principal coordir^ate arialysis of all pairwise Fst values for all 
spatial and temporal sample collections for Siganus spinas. Both axes 
combined explain 66.6% of the total variation. Site codes for the region- 
wide samples are as Table 1 , with Guam recruit samples labeled as GuR. 
Samples are color coded by island. Note that all recruit samples cluster 
close to the Guam/Saipan adult samples. 

STRUCTURE, simulation and empirical evidence has shown 
that the inclusion of sampling information does not lead 
to detection of genetic structure when genetic structure is 
not present (Hubisz et al. 2009). Additionally, the trend of 
isolation-by-distance, which we detected across the region, 
refutes the hypothesis of a single panmictic population be- 
tween the Marshall Islands and the Philippines. 

Our study is the first to document significant genetic struc- 
ture between the southern Mariana Islands and its nearest 
neighbors in Micronesia, and this has potentially profound 
considerations for future management of reef fish popula- 
tions in the region. Previous genetics-based studies within 
the Pacific Ocean basin have generally focused on connectiv- 
ity from a phylo geographic perspective. These studies used 
more highly conserved allozyme and mtDNA markers, which 
reflect historical patterns of connectivity that may not be 
representative of contemporary patterns of gene flow. Con- 
trasting patterns of genetic structure were nevertheless ob- 
served across a range of fish species, although most of the 
cited studies reported an absence of genetic structure (Bay 
et al. 2004; Craig et al. 2007; Home et al. 2008; Gaither et al. 
2010). One study, using allozyme markers, did however find 
strong genetic differentiation across the Pacific, in the sur- 
geonfish Acanthurus triostegus (Planes and Fauvelot 2002). 
Populations from Guam, Philippines, Palau, and the Great 
Barrier Reef (GBR) were all genetically differentiated from 
each other but formed a western Pacific clade more simflar 
to each other than the rest of the Pacific. In the only other ex- 
ample of connectivity assessment across Micronesia, Rhodes 
et al. (2003) found strong genetic structure in the camouflage 
grouper Epinephelus polyphekadion using three micros at ellite 
loci. Cluster analysis identified three regions of genetic differ- 



entiation: Pohnpei and the Marshall Islands; New Caledonia 
and the GBR; and Palau. The separation of Pohnpei and the 
Marshall Islands from Palau is also evident in S. spinus (this 
study), which could indicate a consistent barrier to disper- 
sal between East and West Micronesia. However, inter annual 
variation was observed within Pohnpei, which may suggest 
populations of this species are not temporally stable. 

In rabbitfish, large variations in patterns of population 
structure have been observed between species. Using mtDNA, 
Iwamoto et al. (2009) found no evidence of genetic structure 
in S. spinus among islands in Okinawan waters separated 
by 430 km. In contrast, S. guttatus sampled concurrently 
did exhibit genetic structuring. The authors attribute this 
disparity to the early life-history characteristics of S. spinus 
with their larvae having greater dispersive potential than S. 
guttatus. A similar result was obtained when S. argenteus, 
a species with early life-history characteristics similar to S. 
spinus (Woodland 1990; Iwamoto et al. 2009), was exam- 
ined for genetic structure along the Philippine coast, along 
with S.fuscescens (Magsino and Juinio-Menez 2008). Genetic 
structure was present in S. fuscescens, but not in S. argenteus. 
Again, this was attributed to early life-history characteris- 
tics, and also oceanographic conditions. Further supporting 
this conclusion, a molecular systematic comparison of eight 
species of rabbitfish found that S. argenteus and S. spinus 
exhibited the lowest intraspecific genetic differentiation of 
the genus (Lemer et al. 2007). Our results also demonstrate 
high connectivity over large distances, as demonstrated by 
the delineation of the East Micronesia population that spans 
a distance of more than 2000 km. With a pelagic larval dura- 
tion of 32 days, strong swimming ability and large settlement 
size, S. spinus undoubtedly possess traits that appear to favor 
long-distance dispersal. However, simple life-history traits 
have shown to be poor predictors of both genetic differen- 
tiation and self-recruitment (Bay et al. 2006; Aim any et al. 
2007; Weersing and Toonen 2009; but see Riginos and Vic- 
tor 2001). Also, larval behavior, such as swimming, school- 
ing, and natal homing could actively reduce dispersal (Selkoe 
et al. 2008). While little is known of S. spinus larval behav- 
ior, presettlement stage fish have been observed on Guam 
schooling offshore prior to recruitment onto reefs (Kami and 
Ikehara 1976). 

Oceanographic models have been used to predict levels of 
dispersal and have shown in many cases that prevailing cur- 
rent regimes are a driving mechanism of connectivity, and, 
thus, the genetic structure of marine populations (Glig and 
Hilbish 2003; Galindo et al. 2006; Galarza et al. 2009; White 
et al. 2010). The southern Mariana Islands are located in the 
path of the North Equatorial Current (NEC) , which flows pre- 
dominately northwestward at approximately 0.1-0.3 msec~^ 
(Fig. 1). In contrast, all our other sampling locations (except 
PNG) are located in lower latitudes and are under the influ- 
ence of the eastward flowing Equatorial Countercurrent. A 
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recent biophysical model constructed using satellite-derived 
broad-scale current regimes to predict population connectiv- 
ity of coral larvae, revealed complete isolation of the southern 
Mariana Islands from the rest of the Pacific when the pelagic 
larval duration of simulated larvae was set to 30 days (Treml 
et al. 2008). This pattern was consistent even when accounting 
for oceanographic variability as a result of El Nino-La Nina 
events. Their model also predicted high levels of connectivity 
across the Micronesian islands, which strongly agrees with the 
large-scale spatial structure we report here; suggesting ocean- 
scale current flows are a significant influence on patterns of 
connectivity in S. spinus. At the local scale, transient eddies 
have been observed forming in the lee of Guam, generated by 
the NEC. These eddies maybe sufficiently energetic to return 
larvae to their natal reefs (Wolanski et al. 2003), revealing a 
potential mechanism for enhanced self-recruitment within 
the southern Mariana Islands. A biophysical model recently 
compiled for Guam has also provided evidence in support of 
self-entrainment through island eddy formation (A. Halford, 
pers. comm.). 

In contrast to many similar studies, we also investigated 
the temporal strength of our observed spatial structuring and 
found no genetic differences between new cohorts and adult 
S. spinus on Guam. We also found no evidence of sweepstakes 
reproduction as a mechanism for structuring populations: 
individuals within recruitment events did not (on average) 
share more alleles with other recruits or with adults from 
the same site. This suggests the southern Mariana Islands 
should be managed as a single demo graphically connected 
population. Temporal stability amid low levels of gene flow 
suggests persistent self-recruitment may be responsible for 
maintaining the distinct genetic structure found across the 
southern Mariana Islands (Swearer et al. 2002). Further in- 
direct evidence for self-recruitment comes from the PCoA 
of pairwise Fsr values for all sample groups. When grouped, 
allele frequencies from the recruit samples were not differen- 
tiated from adult samples from the southern Mariana Islands 
but showed significant differences when compared to adult 
populations from both the East Micronesia and West Pacific 
populations. Direct genetic evidence for self- recruitment in 
marine reef fish is steadily increasing. Studies using parentage 
analysis and assignment tests have shown locally produced 
larvae can be a primary source of population replenishment 
(Jones et al. 2005; Planes et al. 2009; Christie et al. 2010). De- 
spite Bayesian assignment tests failing to assign most recruits 
to a single population of origin, the seven recruits that could 
be assigned all originated from the southern Mariana Islands 
population, providing some affirmation of self- recruitment 
as a significant process in structuring the population. The 
low power of assignment is likely a result of the low number 
of markers used and their relatively low variability. Empir- 
ical and simulation studies have shown low polymorphic 
loci provide less accurate assignments than high polymor- 



phic loci (Waples and Gaggiotti 2006; Saenz-Agudelo et al. 
2009). Reduced polymorphism is a common problem with 
loci originally designed for a different species (EUegren et al. 
1997; Neff and Gross 2001), and we understand that these 
limitations may have hindered our ability to assign recruits 
and detect genetic structure using individual-based analyses 
at finer resolutions, both spatially and temporally. Unfortu- 
nately, financial and logistical constraints prevented the de- 
velopment of species-specific markers in this study and it is 
anticipated that further insights may be gained from the use 
of additional loci with greater polymorphism (Kalinowski 
20024). However, the life-history characteristics of S. spinus 
may preclude the use of parentage analysis. Rabbitfish are 
fast growing and subsequently suffer high natural mortality 
(Ntiba and Jaccarini 1988; Grandcourt et al. 2007). Siganus 
spinus can also reach sexual maturity within approximately 
6 months (A. Halford, unpubl. data). Such a dynamic de- 
mographic structure may dramatically reduce the chance of 
finding parent-offspring pairs. Nonetheless, given that, un- 
der certain conditions, very few successful migrants per gen- 
eration can lead to genetic homogeny among populations, 
our results provide strong evidence for treating the southern 
Mariana Islands as an isolated, predominately self- recruiting 
population (Mills and Allendorf 1996). 

The results presented in this study may have significant 
implications for the management of Guam's S. spinus fish- 
ery. As a species subjected to high fishing pressure and with 
historically high variance in recruitment, S. spinus has been 
recognized on Guam as a species in need of conservation. Our 
results suggest a greater emphasis on effective local manage- 
ment of fish stocks is required, because recruitment may be 
directly linked to the standing stock of local adult S. spi- 
nus. Moreover, suitable habitat for S. spinus is limited in the 
Northern Mariana Islands as these islands are much younger 
geologically and do not contain significant reef flat habitat 
(Houk and Starmer 2010), and are unlikely to be acting as a 
larval source for this species. Currently, there are five Marine 
Preserves on Guam, three of which contain reef flat habitat 
suitable for S. spinus. Each of these three preserves contains 
greater S. spinus biomass than adjacent fished sites (A. Hal- 
ford, unpubl. data). However, the preserves are opened sea- 
sonally to allow fishing for newly settled S. spinus recruits. 
The effects of directly fishing recruitment pulses (as in Pauly 
et al. 1998) within the preserves are unknown, but the results 
presented here suggest adequate protection of adult spawn- 
ing stocks may be essential for building resilience in the local 
population of S. spinus and to help ensure the future sustain- 
ability of the fishery. 

Further research into connectivity within the Mariana 
archipelago region is clearly needed. If the results presented 
here are mirrored in other species, then managing for future 
resilience of local reef fish populations will have to explicitly 
consider the strength of the stock-recruitment relationship. 
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An integral part of this relationship is the magnitude of lar- 
val exchange between islands within the southern Mariana 
archipelago (Guam and Saipan). While clearly enough to 
maintain genetic homogeneity, the actual extent of larval ex- 
change between islands remains unknown. Such information 
is critical for providing effective conservation of exploited 
marine resources. 
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Table Al. Characteristics of five microsatellite loci from 1 8 collections of Siganus spinus (>60 mm FL) from the western Pacific. Given are: Site code 
and sample size (A/), number of alleles {A), observed heterozygosity (Ho), expected heterozygosity (H^), and inbreeding coefficient (Fis)^ . 



Sample code (A/) 




Sfus-113 


Sfus-56 


Sfus-8 


Sfus-9 


Sfus-98 


PId 


A 


13 


2 


19 


2 


2 


88 


Ho 


0.824 


0.011 


0.852 


0.264 


0.138 




He 


0.851 


0.011 


0.919 


0.229 


0.185 




F^s 


0.038 


0.000 


0.080 


-0.147 


0.262 


PAb 


A 


15 
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18 


2 


2 


50 


Ho 


0.804 


0.040 


0.786 


0.146 


0.163 




He 


0.858 


0.039 


0.918 


0.170 


0.215 




F^s 


0.074 


-0.010 


0.156 


0.152 


0.250 


PAk 


A 


12 


2 


15 


2 


2 


40 


Ho 


0.917 


0.050 


0.750 


0.179 


0.175 




He 


0.845 


0.049 


0.890 


0.242 


0.160 




F^s 


-0.070 


-0.013 


0.171 


0.271 


-0.083 


YAw 


A 


9 


2 


13 
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2 


26 


Ho 


0.760 


0.154 


0.692 


0.192 


0.269 




He 


0.793 


0.142 


0.894 


0.233 


0.286 




F^s 


0.062 


-0.064 


0.244 


0.194 


0.079 


PNk 


A 


8 


2 


9 


2 


2 


9 


Ho 


0.875 


0.111 


0.667 


0.222 


0.111 




He 


0.820 


0.105 


0.852 


0.198 


0.401 




F^s 


0.000 


0.000 


0.273 


-0.067 


0.750 


GUc 


A 


9 


2 


12 


1 


3 


26 


Ho 


0.846 


0.038 


0.583 


0.000 


0.269 




He 


0.798 


0.038 


0.852 


0.000 


0.240 




F^s 


-0.041 


0.000 


0.334 


NA 


-0.101 



Continued 



3210 



(c) 2012 The Authors. Published by Blackwell Publishing Ltd. 



M. A. Priest et al. 
Table A1. Continued 
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Sample code (A/) 




Sfus-n3 


Sfus-56 


Sfus-8 


Sfus-9 


Sfus-98 






8 


2 


1 2 


1 


2 


39 


''0 


0.744 


0.026 


0.447 


0.000 


0.205 




Hr 


0.730 


0.025 


0.896 


0.000 


0.1 84 




'IS 


-0.005 


0.000 


0.511 


NA 


-0.101 


GU| 


A 


12 


2 


1 5 


2 


3 


67 




0.800 


0.01 5 


0.621 


0.01 5 


0.258 




He 

' 'E 


0.772 


0.01 5 


0.889 


0.01 5 


0.275 




'IS 


-0.028 


0.000 


0.310 


0.000 


0.072 


GUp 


/\ 


1 1 


2 


1 2 


1 


2 


28 


''0 


0.821 


0.107 


0.583 


0.000 


0.250 




Hr 


0.783 


0.101 


0.856 


0.000 


0.21 9 




'^IS 


-0.032 


-0.039 


0.337 


NA 


-0.125 


GUt 


/\ 


1 2 


2 


1 3 


2 


3 


87 


Hr. 

''0 


0.720 


0.048 


0.579 


0.01 1 


0.184 




Hr 


0.769 


0.046 


0.895 


0.01 1 


0.1 87 




F,c 

1 IS 


0.070 


-0.01 8 


0.359 


0.000 


0.020 


SA^ 


A 


8 


2 


1 1 


1 


3 


1 8 


Hr. 


0.750 


0.1 1 1 


0.647 


0.000 


0.333 




Hr 

' 'E 


0.768 


0.105 


0.808 


0.000 


0.285 




'^IS 


0.055 


-0.030 


0.228 


NA 


-0.140 


SA, 




10 


2 


1 3 


1 


3 


56 


Hr. 


0.745 


0.055 


0.648 


0.000 


0.268 




Hr 


0.753 


0.053 


0.891 


0.000 


0.236 




'^IS 


0.01 9 


-0.01 9 


0.281 


NA 


-0.126 






9 


2 


1 0 


2 


2 


19 


Hr. 

''0 


0.684 


0.105 


0.500 


0.053 


0.368 




Hr 


0.71 5 


0.1 00 


0.859 


0.051 


0.361 




1 IS 


0.070 


-0.029 


0.444 


0.000 


0.008 


CHp 




1 2 


1 


1 6 


2 


2 


92 


Hr. 


0.800 


0.000 


0.525 


0.258 


0.228 




Hr 

' 'E 


0.789 


0.000 


0.91 0 


0.273 


0.202 




'^IS 


-0.009 


NA 


0.428 


0.058 


-0.124 


CHx 


/\ 


12 


1 


1 2 


2 


2 


51 


Hr. 

''0 


0.780 


0.000 


0.405 


0.362 


0.21 6 




Hr 


0.799 


0.000 


0.876 


0.347 


0.1 92 




'^IS 


0.034 


NA 


0.547 


-0.032 


-0.1 1 1 






1 2 


1 


1 5 


2 


2 


1 14 


Hr. 


0.738 


0.000 


0.524 


0.167 


0.257 




Hr 


0.786 


0.000 


0.902 


0.181 


0.237 




Fis 


0.066 


NA 


0.423 


0.086 


-0.080 


POn 


A 


13 


1 


14 


2 


2 


108 


Ho 


0.800 


0.000 


0.522 


0.279 


0.241 




H, 


0.799 


0.000 


0.898 


0.254 


0.265 




F,s 


0.004 


NA 


0.424 


-0.095 


0.097 


MAa 


A 


10 


1 


12 


2 


2 


53 


Ho 


0.792 


0.000 


0.413 


0.113 


0.377 




He 


0.816 


0.000 


0.878 


0.107 


0.306 




F,s 


0.038 


NA 


0.537 


-0.051 


-0.224 



^Numbers in boldface type indicate significant departure from Hardy-Weinberg equilibrium at the 0.05 level after Bonferroni correction. 
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Table A2. Characteristics of six microsatellite loci from samples of Siganus spinas recruits and adults (>100 mm FL) from Guam. Given are: Collection 
code and sample size (A/), number of alleles {A), observed heterozygosity (Ho), expected heterozygosity (H^), and inbreeding coefficient (F/sJ^ . 



Sample code (A/) 




Sfus-113 


Sfus-56 


Sfus-8 


Sfus-9 


Sfus-98 


Sfus-21 


AchN7 


A 


10 


1 


1 1 


1 


3 


6 


20 


' '0 


0.737 


0.000 


0.471 


0.000 


0.100 


0.750 




' 'E 


0.706 


0.000 


0.841 


0.000 


0.096 


0.71 3 




'IS 


0.070 


NA 


0.554 


NA 


-0.01 3 


-0.027 


AgaN7 




6 


2 


12 


2 


2 


5 


13 


''0 


0.846 


0.077 


0.538 


0.077 


0.1 54 


0.692 




' 'E 


0.722 


0.074 


0.882 


0.074 


0.142 


0.707 




Fi5 


-0.133 


0.000 


0.423 


0.000 


-0.044 


0.061 


AgaSS 


/\ 


9 


2 


13 


1 


2 


7 


20 


Hr. 


0.900 


0.050 


0.700 


0.000 


0.1 00 


0.700 




He 


0.806 


0.049 


0.858 


0.000 


0.095 


0.748 






-0.091 


0.000 


0.208 


NA 


-0.027 


0.089 


CocA 


A 


9 


2 


12 


2 


3 


7 


24 


Hr. 

' '0 


0.875 


0.042 


0.565 


0.042 


0.333 


0.542 




He 


0.803 


0.041 


0.857 


0.041 


0.288 


0.61 6 




'IS 


-0.069 


0.000 


0.394 


0.000 


-0.136 


0.142 


EAgN7 


A 


8 


2 


8 




2 


6 


16 


Hr. 

' '0 


1 .000 


0.125 


0.429 


0.000 


0.125 


0.81 3 




He 

' 'E 


0.801 


0.1 1 7 


0.829 


0.000 


0.1 1 7 


0.693 




' IS 


-0.218 


-0.035 


0.582 


NA 


-0.035 


-0.140 


EAgS8 


/\ 


10 




1 3 


2 


3 


6 


21 


Hr. 
• '0 


0.762 


0.000 


0.429 


0.048 


0.238 


0.762 




He 

' 'E 


0.787 


0.000 


0.887 


0.046 


0.214 


0.662 




Fk 
' IS 


0.056 


NA 


0.534 


0.000 


-0.087 


-0.127 


GovA 




8 


2 


13 




3 


7 


26 


Hr. 

' '0 


0.769 


0.038 


0.542 


0.000 


0.423 


0.577 




He 


0.737 


0.038 


0.891 


0.000 


0.341 


0.61 6 




Fk 
' IS 


-0.024 


0.000 


0.461 


NA 


-0.222 


0.083 


GnvOR 


/\ 


8 


2 


10 


2 


2 


7 


20 


Hr. 

''0 


0.737 


0.050 


0.444 


0.050 


0.100 


0.550 




He 

' 'E 


0.745 


0.049 


0.847 


0.049 


0.095 


0.626 




F,5 


0.113 


0.000 


0.556 


0.000 


-0.027 


0.147 


GovS8 


/\ 


9 


2 


1 1 


1 


2 


6 


26 


Hr. 
• '0 


0.808 


0.038 


0.478 


0.000 


0.1 1 5 


0.538 




He 

' 'E 


0.752 


0.038 


0.882 


0.000 


0.174 


0.531 




Fk 
' IS 


-0.054 


0.000 


0.541 


NA 


0.353 


0.006 


IpaA 




1 1 


2 


14 


2 


3 


6 


42 


Hr. 

' '0 


0.805 


0.024 


0.61 5 


0.024 


0.167 


0.732 




He 
• 'E 


0.761 


0.024 


0.889 


0.024 


0.194 


0.682 




Fk 
' IS 


-0.006 


0.000 


0.375 


0.000 


0.151 


-0.013 


lpaN7 


/\ 


8 


2 


10 


2 


2 


5 


13 


Hr. 

' '0 


0.727 


0.077 


0.500 


0.077 


0.077 


0.462 




He 


0.802 


0.074 


0.833 


0.074 


0.074 


0.444 




Fis 


0.299 


0.000 


0.489 


0.000 


0.000 


0.000 


IpaSS 


^ 


10 


1 


10 


1 


2 


6 


19 


Ho 


0.789 


0.000 


0.368 


0.000 


0.316 


0.579 




He 


0.755 


0.000 


0.853 


0.000 


0.266 


0.683 




Fis 


-0.019 


NA 


0.586 


NA 


-0.161 


0.178 


PagA 


>A 


9 


2 


6 


1 


2 


5 


11 


Ho 


0.818 


0.091 


0.500 


0.000 


0.364 


0.545 




He 


0.789 


0.087 


0.765 


0.000 


0.298 


0.636 




F,s 


0.011 


0.000 


0.468 


NA 


-0.177 


0.189 
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Sample code (A/) 




Sfus-113 


Sfus-56 


Sfus-8 


5fus-9 


5fus-98 


Sfus-21 


PagSS 




9 


1 


1 3 


1 


3 


5 


22 


''0 


0.864 


0.000 


0.667 


0.000 


0.136 


0.636 




Hr 


0.774 


0.000 


0.872 


0.000 


0.1 29 


0.661 




'IS 


-0.093 


NA 


0.299 


NA 


-0.033 


0.061 


RioJS 


A 


8 




10 




2 


7 


1 5 




0.733 


0.000 


0.333 


0.000 


0.200 


0.600 




He 

' 'E 


0.773 


0.000 


0.880 


0.000 


0.180 


0.562 




'IS 


0.086 


NA 


0.642 


NA 


-0.077 


-0.033 


RioOS 


/\ 


5 


1 


8 






5 


8 




0.800 


0.000 


0.375 


0.000 


0.000 


0.500 




Hr 


0.740 


0.000 


0.836 


0.000 


0.000 


0.609 




'^IS 


0.398 


NA 


0.596 


NA 


NA 


0.243 


RioS8 




1 2 


1 


1 4 


1 


3 


5 


48 


Hr. 


0.771 


0.000 


0.548 


0.000 


0.292 


0.563 




Hr 


0.733 


0.000 


0.894 


0.000 


0.259 


0.641 




1 IS 


-0.041 


NA 


0.478 


NA 


0.463 


0.1 34 


Tan A 




10 


2 


1 3 




3 


7 


33 


Hr. 


0.71 9 


0.030 


0.594 


0.000 


0.1 52 


0.697 




Hr 

' 'E 


0.793 


0.030 


0.900 


0.000 


0.142 


0.685 




1 IS 


0.149 


0.000 


0.377 


NA 


-0.053 


-0.003 


TanJB 


A 


8 


1 


8 




3 


7 


1 9 


H^ 


0.632 


0.000 


0.474 


0.000 


0.263 


0.789 




Hr 


0.673 


0.000 


0.845 


0.000 


0.234 


0.71 6 




'^IS 


0.089 


NA 


0.461 


NA 


-0.098 


-0.076 


TanN7 


A 


8 


1 


7 


1 


2 


4 


9 


Hr. 


0.778 


0.000 


0.556 


0.000 


0.222 


0.667 




Hr 


0.784 


0.000 


0.796 


0.000 


0.1 98 


0.574 




Fis 


0.067 


NA 


0.355 


NA 


-0.067 


-0.103 


TanOS 


A 


8 


1 


11 


1 


2 


6 


23 


Ho 


0.826 


0.000 


0.545 


0.000 


0.000 


0.652 




H, 


0.787 


0.000 


0.870 


0.000 


0.083 


0.680 




F,s 


-0.027 


NA 


0.425 


NA 


1.000 


0.063 


TangSS 


A 


8 


1 


11 


2 


2 


7 


19 


Ho 


0.789 


0.000 


0.526 


0.053 


0.105 


0.737 




He 


0.770 


0.000 


0.875 


0.051 


0.188 


0.652 




F,s 


0.002 


NA 


0.421 


0.000 


-0.116 


-0.103 



^Numbers in bold typeface indicate significant departure from Hardy-Weinberg equilibrium at the 0.05 level after Bonferroni correction. 
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